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INCORPORATION  OF  DRIFT  LOSS  CONE  EFFECTS  ON  THE  TRAPPING  OF  AN  ARTIFICIAL 
RADIATION  BELT  INTO  LANL'S  MODELING  CAPABILITY 

G.S.  Cunningham  and  M.M.  Cowee 
Los  Alamos  National  Laboratory 
Los  Alamos,  NM 

Abstract 

Los  Alamos  National  Laboratory  has  developed  an  end-to-end  model  that  quantifies  the  total  ionizing 
dose  that  a  shielded  component  on  a  space  asset  will  accrue  due  to  interactions  with  trapped  beta- 
decay  electrons  produced  after  a  high-altitude  nuclear  explosion.  The  model  consists  of  three  major 
components:  initial  trapping,  spatio-temporal  evolution  of  the  trapped  population,  and  interaction  of 
the  evolving  trapped  population  with  a  space  asset  to  produce  total  ionizing  dose.  The  initial  trapping 
module  calculates  the  fraction  of  betas  produced  by  radioactive  decay  of  nuclear  fission  products  that 
will  be  trapped  in  the  magnetic  field  surrounding  Earth  for  at  least  one  drift  around  Earth.  Previously,  the 
initial  trapping  calculation  assumed  a  centered  dipole  magnetic  field  model  whereas  we  know  that  Earth's 
magnetic  field  is  offset  from  its  geographic  center  and  has  higher-order  non-dipole  terms.  Both  the  offset 
and  the  higher-order  terms  are  captured  in  the  International  Geomagnetic  Reference  Field  (IGRF)  internal 
field  model.  In  a  centered  dipole,  the  electrons  that  are  trapped  on  a  magnetic  field  line  at  one  longitude 
will  be  trapped  at  all  other  longitudes  as  they  drift  around  Earth  due  to  azimuthal  symmetry  of  the  field, 
but  the  offset  and  higher-order  terms  of  the  IGRF  field  affects  the  drift  of  the  trapped  electrons  and  causes 
some  that  are  trapped  at  one  longitude  to  be  lost  at  other  longitudes.  These  electrons  are  called  "quasi- 
trapped"  and  are  said  to  be  in  the  drift  loss  cone  since  they  will  be  lost  during  their  first  drift  around  Earth. 
In  this  paper  we  report  on  the  effect  of  incorporating  the  drift  of  electrons  around  Earth  using  the  IGRF 
field  on  the  trapping  fraction  as  compared  to  using  a  centered  dipole.  We  find  that  the  new,  more  accurate, 
methodology  produces  higher  trapping  fractions  for  bursts  above  the  South  Atlantic  Anomaly,  a  region 
of  low  magnetic  field  strength,  and  lower  trapping  fractions  for  bursts  in  regions  with  higher  magnetic 
field  strength  than  expected  from  a  centered  dipole. 


Introduction 

The  high  altitude  nuclear  test  shots  of  the  1 950's  and 
60's  demonstrated  that  high  energy  beta-decay  electrons 
(betas),  produced  by  radioactive  decay  of  fission  frag¬ 
ments,  can  become  trapped  in  Earth's  magnetic  field  and 
form  an  artificial  radiation  belt  or  enhance  a  pre-existing 
natural  belt  [1],  Los  Alamos  National  Laboratory  (LANL) 
has  developed  an  end-to-end  model  that  quantifies  the 
total  ionizing  dose  that  a  shielded  component  on  a  space 
asset  will  accrue  due  to  interactions  with  such  radiation 
belts.  The  LANL  end-to-end  model  consists  of  three  major 
components:  initial  trapping,  spatio-temporal  evolution 
of  the  trapped  population,  and  interaction  of  the  evolv¬ 
ing  trapped  population  with  a  shielded  component  on 
a  space  asset  to  produce  total  ionizing  dose.  The  initial 
trapping  is  computed  by  the  Electron  Source  Model 
(ESM).  The  ESM  uses  a  fission  fragment  debris  model 


that  allows  some  of  the  debris  to  travel  along  field  lines, 
emitting  electrons  isotropically  along  the  way.  Electrons 
that  have  velocity  vectors  that  are  nearly  parallel  to  the 
local  magnetic  field  line  direction  will  travel  down  the 
field  line  to  the  near-surface  atmosphere  and  be  lost, 
whereas  electrons  that  have  velocity  vectors  nearly 
perpendicular  to  the  field  direction  will  "mirror"  in  the 
field  and  be  trapped.  The  gradient  and  curvature  of  the 
magnetic  field  cause  the  electrons  to  drift  around  Earth. 
We  have  recently  improved  the  ESM  so  that  it  can  make 
use  of  the  International  Geomagnetic  Reference  Field 
(IGRF),  which  includes  the  tilt  and  offset  of  the  dipole 
portion  of  the  field  as  well  as  the  higher-order  terms  that 
can  be  modeled  as  a  sum  of  spherical  harmonics  whose 
coefficients  vary  slowly  with  time,  whereas  the  ESM  previ¬ 
ously  assumed  that  the  surrounding  magnetic  field  was  a 
centered  dipole  [2].  Furthermore,  we  also  now  compute  a 
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complete  set  of  closed  drift-shells,  i.e.,  drift  trajectories  for 
which  the  electron's  minimum  altitude  is  above  1 20  km  at 
all  longitudes  and  can  therefore  assume  to  be  trapped  for 
at  least  one  drift  around  Earth  (longitudinal  drift  is  caused 
by  the  gradient  and  curvature  of  the  magnetic  field  and 
takes  on  the  order  of  1 0's  of  minutes  for  an  MeV  electron). 
Combining  the  two  new  capabilities  allows  us  to  compute 
the  effect  of  drifting  in  an  IGRF  field  on  the  trapping  frac¬ 
tion,  which  we  report  on  here. 

First,  we  review  the  basic  particle  motions  in  a  mag¬ 
netic  field  and  their  associated  invariants.  Next  we  discuss 
the  new  capability  that  has  been  introduced  into  the  ESM, 
the  calculation  of  the  drift  shells,  and  the  mapping  of  the 
electrons  produced  by  the  ESM  onto  the  set  of  closed 
drift  shells.  Finally,  we  compare  the  trapping  fractions 
obtained  under  various  assumptions  and  draw  some 
conclusions. 

Background 

Charged  particles  with  a  velocity  v  in  a  magnetic 
field  B  undergo  three  periodic  motions:  gyration  about 
a  magnetic  field  line  due  to  the  Lorentz  qvxB  force, 
movement  along  a  magnetic  field  line  due  to  the  compo¬ 
nent  of  velocity  that  is  parallel  to  the  magnetic  field  vector 
(with  "mirroring"  at  some  altitude  due  to  conservation  of 
the  1st  invariant,  discussed  in  the  next  paragraph),  and 
drift  around  Earth  due  to  the  gradient  and  curvature  of 
the  magnetic  field.  Associated  with  each  periodic  motion 
is  an  adiabatic  invariant  that  is  conserved  so  long  as  there 
are  no  collisions  and  variations  in  the  electric  or  magnetic 
field  occur  on  time  scales  that  are  long  compared  to  the 
periodic  motion. 

The  gyro  motion  is  associated  with  the  first  adiabatic 
invariant  p  =  p~  sin2(or)/(2m0fi),  where  p  is  the  to¬ 
tal  momentum,  a  is  the  local  pitch-angle,  i.e.,  the  angle 
of  the  velocity  vector  with  respect  to  the  magnetic  field 
direction,  m0  is  mass,  and  B  is  magnetic  field  intensity. 
Conservation  of  p  as  a  particle  moves  along  a  field  line 
toward  the  mirror  point  implies  that  a  increases  until 
all  of  the  momentum  is  perpendicular  to  the  field  line, 
i.e.,  a  =90°,  at  which  point  the  particle  begins  to  travel 
in  the  opposite  direction.  Since  a  varies  along  the  field 
line,  the  value  of  a  at  the  point  where  the  intensity  of 
the  magneticfield  is  minimized  is  often  used  and  is  called 


the  equatorial  pitch-angle,  a  ,  since  the  magnetic  field 
intensity  is  minimized  along  any  field  line  at  the  magnetic 
equator  for  a  dipole  field. 

The  bounce  motion  is  associated  with  the  second 
adiabatic  invariant,  J  =  2j  pcos(a)ds,  where  sm  is 
the  arclength  associated  with  the  southern  hemisphere 
mirror  point  and  nm  is  the  arclength  associated  with  the 
northern  hemisphere  mirror  point.  We  use  a  related  ver¬ 
sion  of  the  invariant  K  =  ^Bm  I  =  J  -\j  B,„  ~  B(s)ds, 
where  B(sm)  =  B(nm)  =  Ba  is  the  mirror-point  magnetic 
field  intensityand  isalsoan  invariant  of  motion  underthe 
assumption  that  the  total  momentum  doesn't  change. 

As  the  particle  drifts  longitudinally  around  Earth,  it 
drifts  onto  field  lines  that  preserve  the  first  two  invariants, 
or  equivalently,  Bm  and  K.  At  each  longitude  there  is  a 
unique  field  line  with  one  minimum  value  that  has  the 
particle's  Bm  and  K.  The  set  of  field  lines  that  have  this 
property  are  called  the  drift-shell  of  the  particle.  If  we 
trace  the  field  lines  to  120  km  altitude  in  the  northern 
hemisphere,  the  set  of  points  defines  a  contour  n  inside 
of  which  we  can  integrate  the  magnetic  flux  that  points 
outward,  ®  =  |  B  dS ■  The  quantity  ®  is  an  invari- 
ant  of  motion  associated  with  the  drift  motion  of  the 
particle.  A  related  quantity  that  we  use  is  the  Roederer 
L*  =  2ttBeR2e  /®  [3].  If  the  magnetic  field  is  a  centered 
dipole,  then  r(0)  =  L  sirr(6*)  defines  the  field  lines 
and  L*  —  L,  where  9  is  the  polar  angle  from  the  north 
magnetic  pole,  also  called  the  colatitude.  The  approach 
that  we  use  to  compute  the  drift-shells,  also  used  for  the 
results  presented  in  this  paper,  is  detailed  in  [4]. 

In  a  centered  dipole,  the  value  of  a  for  a  particle 
that  mirrors  at  the  120  km  atmosphere  cutoff  altitude 
is  called  the  bounce  loss  cone  angle  ablc.  Particles  that 
have  a  value  of  ae  <  ablc  will  mirror  below  1 20  km  and 
are  likely  to  be  scattered  and  lost  through  collisions.  We 
use  120  km  as  a  simplification;  in  reality,  the  altitude  at 
which  particles  are  scattered  and  lost  to  the  atmosphere 
is  energy-dependent,  with  higher  energy  particles  able  to 
reach  lower  altitudes.  We  note  that  for  L  <  1 .0 1 88,  the 
particle's  entire  orbital  path  is  below  1 20  km  altitude  and 
so  ablc  =  90  °. 
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In  an  eccentric  dipole,  i.e.,  one  that  is  offset  and  tilted 
relative  to  Earth's  rotation  axis,  the  dipole  field  lines  are 
azimuthally  symmetric  about  the  axis  of  the  dipole,  but 
are  shifted  and  rotated  relative  to  Earth's  rotation  axis. 
Thus,  the  value  of  ablc  depends  on  longitude.  Because 
aeq  is  a  conserved  quantity  in  a  dipole,  the  maximum 
value  of  ahlc  over  all  longitudes,  adlc  =  max  {ablc  (#?)}, 

<p 

called  the  drift  loss  cone,  is  more  relevant  to  the  trapping 
fraction  since  it  defines  the  set  of  equatorial  pitch-angles 
that  will  remain  trapped  over  an  entire  drift  orbit  around 
Earth.  One  can  analytically  calculate  adlc  as  a  function  of 
L  and  compare  it  to  ablc  for  the  centered  dipole  (Fig.  1). 
Similarly,  one  can  compute  adIc  for  the  IGRF  magnetic 
field  model  (circa  1 960),  although  it  is  important  to  note 
that  aeq  is  not  a  conserved  quantity  in  a  non-dipole  field 
and  therefore  the  aeq  that  is  plotted  in  Fig.  1  for  the  IGRF 
model  is  technically  only  valid  at  a  single  longitude.  A 
more  precise  statement  for  the  IGRF  model  would  be  to 
plot  the  maximum  value  of  the  second  invariant,  K,  for 
which  a  closed  drift-shell  exists;  however,  we  chose  not 
to  do  this  because  the  dependence  of  aeq  on  longitude 
is  relatively  weak  in  Fig.  1  (less  than  a  degree). 

Methodology 

Note  that  if  one  were  to  assume  that  particles  start  to 
interact  strongly  with  the  atmosphere  at  500  km  altitude 
and  thus  become  lost  if  they  continue  to  travel  along  the 
field  line  below  that  point,  then  the  plot  of  ablc  vs  L  for 
a  centered  dipole  in  Fig.  1  (the  red  line)  would  be  shifted 
to  the  right  and  look  more  like  the  plots  of  adlc  using  the 


eccentric  dipole  (blue  line)  and  IGRF  1960  (black  line). 
For  this  reason,  because  we  were  constrained  to  use  a 
centered  dipole  field  model  in  our  previous  work,  we  used 
a  500  km  cutoff  for  the  mirror  point  of  trapped  electrons 
to  mimic  the  drift  loss  cone  effect  when  calculating  the 
trapping  fraction  using  the  ESM.  The  present  work  com¬ 
pares  this  approach  for  calculating  trapping  fractions  to 
two  other  approaches. 

The  first  approach  uses  the  value  of  ablc  for  the  field 
lines  on  which  the  beta-decay  electrons  are  initially 
produced.  The  second  approach  maps  the  bounce- 
trapped  electrons  onto  drift-shells  and  computes  the 
trapping  fraction  as  the  number  of  electrons  that  are 
trapped  for  at  least  one  orbit  around  Earth  divided  by 
the  total  number  of  electrons  produced.  In  previous 
work  [4],  we  detail  how  we  construct  a  comprehensive 
set  of  drift-shells  that  consist  of  field  lines  with  a  single 
minimum  value.  There  may  be  closed  drift-shells  that 
have  field  lines  with  more  than  one  minimum  value, 
called  Shabansky-type  orbits  [5],  but  we  do  not  con¬ 
sider  these  orbits  since  they  only  occur  at  larger  L  and 
it  is  not  clear  that  electrons  on  these  orbits  would  stay 
trapped  for  very  long  [6].  Each  drift-shell  is  associated 
witha  pairofsecondand third  invariants,  {Ki,L*l ),and 
consists  of  a  set  of  field  lines  sampling  a  discrete  set 
of  azimuths.  The  minimum  value  of  the  field  intensity 
along  each  field  line  is  associated  with  a  radius  and 
azimuth  pair,  [rjjk ,  (pijk  ],  at  which  point  we  compute  the 
equatorial  pitch-angle,  a^k. 


—  Bounce  loss  cone 

—  Drift  loss  cone  eccentric  dipole 

—  Drift  loss  cone  IGRF  1960 


Figure  I.The  bounce  loss  cone  angle  (Xblc  for  the  centered  dipole  (red),  drift  loss  cone  angle  CCdlc  for  the  eccentric  dipole  (blue) 
and  drift  loss  cone  angle  (Xdlc  for  the  IGRF  1 960  magnetic  field  model  (black)  as  a  function  of  L. 
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Burst  at  -35°N,  -15°E,  300  km 
IGRF  Magnetic  Field 


Burst  at  -35°N,  -15°E,  300  km 
Eccentric  Dipole  Magnetic  Field 


Y 

Figure  2.  Total  number  of  trapped  electrons  versus  L  using  the  eccentric  dipole  (right)  and  IGRF  (left)  magnetic  field  models, 
for  a  burst  above  the  South  Atlantic  Anomaly  at  300  km  altitude.  Trapping  fractions  were  computed  using  a)  the  local  bounce 
loss  cone  with  1 20  km  atmosphere  (x),  b)  the  local  bounce  loss  cone  with  500  km  atmosphere  (circle),  and  c)  by  identifying  those 
electrons  that  are  on  closed  drift  shells  (triangles). 


The  ESM  outputs  the  number  of  electrons  in  each 
of  the  equatorial  pitch-angle  bins  for  each  of  the  field 
lines  along  which  the  fission  fragment  debris  streams 
and  produces  isotropically  distributed  electrons  at  the 
point  of  production.  The  electrons  are  assumed  to  have 
a  Carter-Reines  distribution  in  energy  [7].  Thus,  the  ESM 
output,  integrated  over  energy,  can  be  represented  as 

N\aeq  r  m  I-  where  m  indexes  over  pitch-angle  and 

\  m  ?  n^rnj 

n  indexes  over  field  line.  Note  that  N  is  not  the  number 
density,  but  the  number  density  could  be  obtained  by 
dividing  N  by  an  appropriate  coordinate  volume  oc¬ 
cupied  by  the  particles  in  a  given  drift-shell.  For  each 
{&m  ’^n’tPn)  we  first  finc*  a  Pa'r  bracketing  closed 
drift-shells  such  that  (rn ,  <pn )  is  in  the  domain  of  a  bilinear 
map  of  the  unit  square  onto  the  four  points  that  come 
from  field  lines  at  neighboring  azimuths  from  each  drift 
shell  and  neighboring  indices  into  L  at  fixed  azimuths, 

Le"  (rij,k’<Pij,k)'  (fijMi’PiJMl)'  (rij+iJc><Pi,j+ u)'and 

(7i  j+ 1 k+vVi  j+ 1 k+ i)-  Next'  we  use  t*ie  same  bilinear  map 
to  interpolate  the  a^k  at  the  four  corners  to  provide  an 
estimate  of  aeq  at  each  of  two  neighboring  indices  into 
K  to  ensure  that  a ,q  is  bracketed  by  these  values.  If  these 
two  types  of  brackets  exist,  then  we  say  that  the  elec¬ 
trons  N(a^ ,rn ,(pn)  have  been  mapped  to  a  drift  shell 
and  are  thus  included  in  the  population  that  is  called 
"permanently  trapped,"  meaning  that  the  electrons 


can  drift  around  Earth  without  hitting  the  atmosphere 
at  1 20  km. 

We  used  the  ESM  to  produce  n(cc *q ,  rn ,  (pn )  for  two 
field  models,  eccentric  dipole  and  IGRF  1 960,  for  a  variety 
of  latitude,  longitudes  and  altitudes  (Table  1 )  intended  to 
quantify  the  importance  of  including  the  drift  loss  cone 
in  calculating  trapping  fractions.  The  IGRF  1960  internal 
field  model  was  added  to  theTsyganenkoT89  [8]  external 
field  model,  which  should  not  have  much  affect  at  low  L. 
For  each  ESM  run,  we  computed  three  trapping  fractions. 
First,  we  summed  N{ae^ ,  r„ ,  (pn )  overall  bins  such  that 
aeq  >  ablc  is  satisfied  for  each  sub-population  of  elec¬ 
trons  on  a  given  field  line,  and  divided  this  number  by 
Ntotal,  the  total  number  of  electrons  that  are  produced 

Table  1 .  ESM  scenarios  and  trapping  fractions  (percent)  for  the 
three  approaches. 


Offset  dipole  IGRF  1960 


Lat 

Lon 

Alt 

Bounce 

500  km 

DLC 

Bounce 

500  km 

DLC 

0 

0 

400 

46 

24 

17 

40 

12 

2 

0 

0 

800 

66 

53 

50 

57 

41 

27 

0 

90 

400 

51 

24 

4 

0 

90 

800 

67 

53 

22 

66 

50 

9 

-35 

-15 

300 

52 

0 

23 

51 

0 

24 

-35 

-15 

500 

57 

22 

23 

55 

22 

23 

-35 

-15 

700 

61 

40 

22 

58 

40 

22 

-35 

-15 

1000 

68 

56 

24 

64 

54 

24 

35 

165 

300 

53 

33 

16 

52 

31 

13 

35 

165 

500 

59 

37 

16 

59 

39 

14 

35 

165 

1000 

70 

58 

21 

68 

57 

18 
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through  radioactive  decay  of  fission  fragment  debris. 
This  trapping  fraction  is  listed  in  the  column  labeled 
"Bounce."  Second,  we  summed  N{a^ ,r„,<Pn )  over  all 
bins  such  that  >  ab^  is  satisfied  for  the  field  line  on 

which  the  electrons  are  tied,  where  a™  is  the  bounce 
loss  cone  assuming  that  the  atmosphere  is  at  500  km 
instead  of  120  km,  and  divided  this  sum  by  N total  t0 
get  the  trapping  fractions  in  the  column  labeled  500 
km  in  Table  1.  Third,  we  attempted  to  find  bracketing 
drift-shells  as  described  above  for  each  bin  (m,  n)  and 
summed  N (a ^  )  for  those  ( m ,  n)  where  we  suc¬ 

cessfully  found  a  bracket  and  can  thus  assign  the  elec¬ 
trons  in  this  sub-population  to  a  closed  drift-shell.  We 
then  divided  this  sum  by  Ntotal  to  get  a  third  trapping 
fraction,  labeled  DLC  in  Table  1. 

Results 

Proper  inclusion  of  the  drift  loss  cone  is  important 
when  the  burst  is  above  the  South  Atlantic  Anomaly 
(SAA),  where  the  magnetic  field  is  especially  weak, 
as  illustrated  in  Fig.  2.  In  this  case,  because  the  burst 
location  is  near  the  equator  and  at  low  altitude  (300 
km),  every  field  line  containing  electrons  in  this  case 
is  entirely  below  500  km,  and  so  the  trapping  fraction 
computed  using  a  bounce  loss  cone  with  altitude=500 
km,  as  we  did  previously,  is  0%.  However,  because  the 
magnetic  field  intensity  is  weak  in  this  area,  the  parti¬ 


cles  produced  on  these  field  lines  will  move  away  from 
Earth's  surface  as  they  drift  and  will  be  associated  with 
closed  drift  shells.  The  trapping  fraction  computed  by 
identifying  which  electrons  are  on  closed  drift  shells 
is  24%,  with  little  difference  between  results  obtained 
using  an  eccentric  dipole  field  model  versus  using  IGRF 
1960. 

At  other  longitudes,  as  in  Fig.  3,  we  find  that  the 
method  for  calculating  trapping  fraction  based  on  using 
a  bounce  loss  cone  with  the  atmosphere  at  500  km  alti¬ 
tude  grossly  over-estimates  the  trapping  fraction.  In  this 
case,  the  high  altitude  (800  km)  means  that  the  bounce 
path  of  many  of  the  electrons  will  be  entirely  above  500 
km  and  counted  as  trapped  using  this  model.  However, 
for  burst  longitudes  where  the  magnetic  field  is  strong, 
the  particles  will  move  closer  to  Earth  as  they  drift,  and 
those  with  smaller  values  of  equatorial  pitch-angle  will 
be  lost  as  they  drift.  Calculating  the  trapping  fraction 
by  identifying  which  particles  are  on  closed  drift  shells 
produces  a  much  smaller  trapping  fraction  in  this  case 
than  the  approach  using  the  500k  altitude  atmosphere 
(the  opposite  of  what  we  discussed  in  the  previous 
paragraph).  Furthermore,  in  this  case,  the  higher-order 
moments  represented  by  IGRF  also  affect  the  trapping 
fraction  since  the  trapping  fraction  computed  using 
IGRF  1 960  is  2.5x  smaller  than  that  computed  using  an 
eccentric  dipole. 
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Figure  3.  Total  number  of  trapped  electrons  versus  L  using  the  eccentric  dipole  (right)  and  IGRF  (left)  magnetic  field  models, 
for  a  burst  far  from  the  South  Atlantic  Anomaly  at  800  km  altitude.  Trapping  fractions  were  computed  using  a)  the  local  bounce 
loss  cone  with  1 20  km  atmosphere  (x),  b)  the  local  bounce  loss  cone  with  500  km  atmosphere  (circle),  and  c)  by  identifying  those 
electrons  that  are  on  closed  drift  shells  (triangles). 
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Conclusions 

We  have  extended  the  Electron  Source  Model  (ESM)  so 
that  it  transports  the  fission  fragment  debris  along  a  finite 
number  of  field  lines  intersected  by  the  burst  volume  in 
the  IGRF  magnetic  field  model,  whereas  previously  a  less 
accurate  centered  dipole  model  was  used.  This  initial 
population  of  electrons  is  transported  to  the  point  along 
the  field  line  where  the  minimum  value  of  magnetic  field 
intensity  is  obtained  and  the  number  of  electrons  as  a 
function  of  equatorial  pitch-angle  is  tallied  for  each  field 
line.  We  further  computed  the  set  of  closed  drift-shells  in 
IGRF  so  that  we  could  determine  which  electrons  would 
be  "permanently  trapped,"  i.e.,  could  drift  entirely  around 
Earth  without  encountering  the  atmosphere  at  120  km. 
This  population  of  particles  is  only  a  fraction  of  those  par¬ 
ticle  whose  bounce  motion  allows  them  to  stay  above  1 20 
km  altitude  on  the  particularfield  line  on  which  they  were 
produced  through  radioactive  decay,  i.e.,  those  particles 
that  are  "bounce-trapped." 

The  trapping  fraction  was  computed  using  three  differ¬ 
ent  methods,  which  revealed  that  there  can  be  substantial 
sensitivity  of  the  trapping  fraction  to  the  latitude,  longitude, 
altitude  and  field  model  that  is  used. The  new  methodology 
that  computes  permanently-trapped  electrons  using  the 
IGRF  produces  larger  trapping  fractions  than  what  is  ob¬ 
tained  using  a  centered  dipole  model  with  a  500  km  atmo¬ 
sphere  for  burst  longitudes  that  have  comparatively  small 
magnetic  field  strength,  as  in  the  South  Atlantic  Anomaly, 
whereas  the  new  methodology  produces  lower  trapping 
fractions  for  bursts  at  longitudes  with  strong  magnetic  fields. 
We  conclude  that  proper  consideration  of  the  drift  loss  cone 
is  important  in  computing  the  trapping  fraction  of  betas 
produced  by  a  high-altitude  nuclear  explosion. 

In  future  work,  we  will  use  the  mapping  of  the  elec¬ 
trons  onto  closed  drift-shells  to  produce  an  initial  condi¬ 
tion  phase  space  density  that  can  be  evolved  forward  in 
time  in  response  to  interactions  with  naturally-occurring 
electromagnetic  waves  as  modeled  by  the  Dynamic  Ra¬ 
diation  Environment  Assimilation  Model  3D  (DREAM3D), 
a  Fokker-Planck  diffusion  code  that  we  have  used  to 
model  the  natural  radiation  belts  [9]. 
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